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ABSTRACT 

This note documents the methods used by the smartphonc apphcation, "GalMass," which has been released on the 
Android Market. GalMass estimates the halo virial mass (Mvir), stellar mass (Mstar), gas mass (Afgas), and galaxy gas 
fraction of a central galaxy as a function of redshift (z < 2), with any one of the above masses as an input parameter. 
In order to convert between Mvir and Mstar (in either direct i on), GalMass uses fi t ting f unc tions that appro ximate 
the abundance matching models of either IConrov fc Wechslei] (pOOl . iMoster et all (pOlOt ). or iBehroozi etall ^VS). 
G alMass uses a a semi-e mpirical fit to observed galaxy gas fractions to convert between Mgtar and Mgas, as outlined 
in iStewart et all (|2009b[) . 



1. INTRODUCTION 

In the LCDM model of structure formation, dark mat- 
ter (DM) halos form by the hierarchical accretion of 
smaller dark matter halos, as well as mor e diffuse accre- 
tion of matter that is not virialized (e.g. IStewart et al.l 
120081 ) . In this model, galaxies reside at the center of dark 
matter halos, such that the mass scale of a galaxy's en- 
vironment can be characterized by the virial mass of the 
dark matter halo in which it resides. This mass scale 
influences a great number of galaxy properties, including 
luminosity, morphology, star formation rate, gas content, 
and color. As such, theoretical formalisms often catego- 
rizes dark matter halo properties as a function of virial 
mass. 

Unfortunately, halo virial masses are difficult to deter- 
mine, empirically. Instead, observations are more likely 
to categorize galaxies by their stellar mass (or luminos- 
ity). As a result, it is often difficult to compare theoreti- 
cal expectations with observations because the mapping 
between stellar mass a nd virial mass is non -trivial, and 
evolves over time (e.g.. IBehroozi et al.ll201Cil ). For exam- 
ple, theoretical measures of the dark matter halo merger 
rate typically describes mergers in terms of the ratio of 
dark matter masses prior to the merger, while observa- 
tional studies often refer to the stellar mass ratio, even 
though a 1 : 3 ratio in dark matter can easily correspond 
to a 1 : 20 or 1 : 2 ratio in s tellar mass, d epending on 
the mass regime in question ()Stewartll2009D . Even in a 
less formal setting than scientific articles, theorists and 
observers often "think" in terms of these different quan- 
tities, making comprehensive discourse of galaxy evolu- 
tion difficult. The goal of GalMass is to convert between 
different galaxy masses quickly and accurately from the 
convenience of a smartphone application. 

A secondary goal of GalMass is to encourage the scien- 
tific community to consider the development of scientific- 
level applications in the growing smartphone market- 
place (e.g., the "Android Market" or the Apple's "App 
Store" ) . While a number of astronomy applications exist 
in these markets, they are typically intended for a public 
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audience. There are very few applications that are devel- 
oped from within the scientific community for productive 
use in the scientific community (examples of such scien- 
tific applications developed by professional astronomers 
include GravLensHD, CosmoCalc, and Exoplanet j3- 

2. ABUNDANCE MATCHING MODELS 

One method for mapping between the dark matter 
virial mass (Mvir) and the corresponding stellar mass 
(Mstar) presumes a monotonic relationship between halo 
mass and stellar mass. That is, so long as a galaxy is 
the central galaxy to its dark matter halo (not a satel- 
lite galaxy), more massive galaxies are presumed to re- 
side in more massive DM halos. In this model, one can 
achieve a mapping between halo mass and stellar mass 
by setting the abundance of halos more massive than 
Afvir, nfi{> Mdm), equal to the abundance of galax- 
ies more massive than Mgtar, ng{> Mgtar)- (Note that 
in order to correctly account for satellite galaxies us- 
ing this method, one must track each satellite back to 
when it was last isolated, matching the abundance of 
its DM halo at that time, instead of using its current 
subhalo mass.) This technique, often referred to as 
"abundance matching" has been shown to reproduce a 
number of galaxy clustering statistics, as well as match- 
ing observed clo se pair counts of galaxi es as a function 
of redshift (e. g.. IConrov et al.ll200 6: Bcr rier et al]|200l 



iShankar et al 



I.. lUonr 
120061: IS 



tewart et al.ll2009aD . 



While a monotonic relationship between halo mass 
and stellar mass cannot be true for every galaxy on 
a case-by-case basis, especially at lower masses where 
galaxy formation is expecte d to be more stochastic (e.g., 
iBovlan-Kolchin et al.l [20111) . in a statistically averaged 
sense it provides a good characterization of the rela- 
tionship between halo mass and stellar mass that must 
be true in order for LCDM halos to reproduce the ob- 
served stellar mass function. Indeed, even when extrap- 
olating abundance matching relations to very low halo 
masses where it is not well-tested, it has been shown to 
produce a plausible (if not robust) model for the dwarf 
:alaxy satellite popu lation around Milky Way size halos 
Bullock et"all[2010l ). 

In order for GalMass to quickly and accurately es- 
timate virial and stellar masses using this abundance 

^ These applications arc available on the iPhone. 



2 



matching formalism, it utilizes analytic fitting functions 
as approximations to published relations. The curve fit- 
ting routine used in constructing these fits is the pub- 
licly available "pymodclfit" routincQ, which minimizes 
in order to determine the best-fit parameters of a 
custom paramctrization to a set of data. Currently, 
GalMass allows for 3 choic es for the abundance match - 
ing model u s ed: t hat o f iConrov &: Wechsleii ()2009[ ). 



iMoster et ahl ()2010[ ). and iBehroozi et al.l (|2010[) . here- 
after CW09, MIO, and BIO respectively. 



2.1. Fitting Functions: \Conrov & Wechslei\l200^ 

As IConrov fc Wechsleii (j2009( ) includes no analytic 
functional forms to is abundance matching results, I 
choose the following functional form to characterize 

M,tar(Mvir): 



(/3-a) 



(1) 



where a and /3 are the low and high mass slopes, and the 
transition between these power laws is described by a 
characteristic halo mass, M2 (with a normalization cor- 
responding to a characteristic stellar mass Mi). Mini- 
mizing to determine the best-fit parameters, I find 
this functional form to be a good approximation for the 
results of CWOfO- Once a fit has been found at each 
epoch (z < 2), I find that the redshift dependence of 
each parameter (or the log of parameters Mi and M2) is 
well-characterized by either a linear, quadratic, or power 
law relation in z, as seen in the right column of Table [TJ 
Because Equation [T] is not easily reversible, I adopt 
the same functional form in order to independently de- 
termine an analytic fit for A/vir(A/star): 



M. 



M2 + Mst, 



(2) 



with Ml now representing a characteristic stellar mass, 
and M2 giving the normalization (halo mass). Again, I 
find this form to provide a good fit to the CW09 data at 
each epoch, with the redshift evolution of each parameter 
given by a quadratic in z. These fits are shown in the 
left column of Table [1] 



2.2. Fitting Functions: iMoster et al]\201d 

For MIO, a fitting function for Mstar(Mvir) is already 
provided, including redshift evolution. For the sake of 
completeness, I include this fitting form (their Equation 
2) here: 



M., 



2Mv 



M™ 
Ml 



/M, 



V Ml 



(3) 

where (m/Afvir)o provides the normalization at the char- 
acteristic halo mass Mi and (3 and 7 control the behavior 

* This python package, created and maintained by Erik 
ToUerud, was also used to creat e the best-fit "fundamen- 
tal curve" of IToUerud et al.l I I2011I) . It is freely available at 
http: / /packages .python. org/PyModelFlt/lndex. html 

^ See i|2.4l for discussion on the goodness of fit of the equations 
presented here. 



at the low and high mass ends. The redshift dependence 
of their free parameters are given in MIO as power laws 
(or linear relations) in l-\-z. These equations (from their 
Equations 23-26 and Table 7) as well as the equation 
above arc included in the right column of Table [TJ 

Again, this fitting form is not easily reversible. In or- 
der to obtain an analytic fit for Mvir(Mstar), I use the 
above fit from MIO to construct a discrete set of data 
points with which to fit Equation[2] This functional form 
provides a good paramctrization to MlO's results, with 
the redshift evolution of each parameter given by the left 
column of Table [TJ 



2.3. Fitting Functions: \Behroozi et al\\201(i 

BIO gives an analytic fit for Afvir(Afstar), but not vice 
versa (opposite from that of MIO, which provides a fit for 
Mstar(Mvir)). For simplicity's sake, I adopt their fitting 
function for the case with no presumed systematic errors 
(their /i = fc = case). Their fitting form for Afvir(A/star) 
(their Equation 21) is as follows: 

logio(Afvir) = logio(Mi) -h/31ogio(M,tar/A^o) , ^ 
(M,t„/Mo)'^ _ 1 (4) 

^l + (M,ta,r/Mo)~T 2 

where Afi and Mq are the normalization and characteris- 
tic mass, /3 controls the low mass end slope, and 5 and 7 
both control the high mass slope. The redshift evolution 
of the various parameters are given by linear relations in 
(1— a), where a is the scale factor, defined by a = l/{l+z) 
such that a value of a = 1 is present day. The various 
redshift evolutions of these parameters (from their Equa- 
tion 22 and Table 2) are all included in the left column 
of Table [U 

As with MIO, this fitting form is not trivially re- 
versible, and I use their fit to construct a discrete data 
set, from which I fit Equation [TJ in order to parameterize 
Mstar(Mvir). Thc fitting function and redshift evolution 
of each parameter is given in the right column of Table 

m 

2.4. Abundance Matching Discussion 

Because in these models I have independently fit mod- 
els for Mvir(Mstar) and Mstar(Mvir) and there is some 
margin of error to each fitting form, I caution that these 
equations are not entirely self-consistent. One may use 
thc above fitting functions to convert from halo mass to 
stellar mass back to halo mass and not return the precise 
Dcgmnmg valuc[3. 

In general, I find that converting from halo mass to 
stellar mass using the above fitting functions tends to 
be more accurate to the original data of CW09. (Thus, 
whenever GalMass shows a small discrepancy between 
derived values for the CW09 model, priority should 
be given to Mstar(Mvir), when the input parameter is 
halo mass). Similarly, for MIO the original function for 
A/star (Mvir) wiU be more accurate, by definition. Since 
thc original fitting function from BIO goes in the oppo- 

^ I originally designed GalMass to recursively call Afstar{Afvir) 
when converting from stellar mass to halo mass, until it produced 
self-consistent results with Mvir (Mstar ) • For reasons that arc un- 
clear, this approach was found to significantly slow down the oper- 
ation, with lag times of 3 — 5 seconds on a typical Android phone. 
In the end, I opted for speed over precision. 



3 



Table 1 

FITTING FUNCTIONS USED BY GALMASS 



Conrov & Wechsler f2009) 


■'■^-'■vir V-'-'-'-star ) 




Fitting Function 






iXTrM'TTi nil "zntir^Ti 

1 N Wi lilCLllZicXUlWll 


Incr, „CMi "1 — n 'VSly'^ — f) 28^1? -1- 1 2 '^^Q 


logio(Afi) = 0.0562^ + .0682 + 9.5 


Characteristic Mass 


logio(A/2) = 0.10222 - 0.0892 + 10.740 


logio(Af2) = 0.3202^ + .0182 + 11.2 


Low Mass Slope 


a. = -0.0183z2 + 0.017262 + 0.3238 


a = .021z4-86 + 3.39 


High Mass Slope 


/3 = 0.12142^ ~ 0.90532 + 3.170 


/3 = .0852 + .36 


Moster et al. f2010") 


Mvir(Mstar) * 


Mstar(Mvir) ^ 


Fitting Function 




Afstar = 2Afvir (iT^)^, 




-1 


Normalization 


logj^(,(A/i) = 0.0222^ - 0.01852 + 12.544 


(m/A/vir)o = 0.0282(1 + 


2)-0-72 




Characteristic Mass 


logio(A/2) = 0.03212^ - 0.1382 + 10.846 


logio(A/i) = 11.884(1 + 


2)0.019 




T .nur A/fnc;c; '^lrtT^e> 

J^U W lV±(XOO OIU L/C 


n- — 473z — 0330 


13 = 0.172 + 1.06 


High Mass Slope 


/3 = -0.58820-489 + 2.64 


7 = 0.556(1 + 2)-0-26 


Behroozi et al. (2010") 


M^i^CM.t.,.') t 


Mstar(Mvir) * 


Fitting Function 


logio(A/vir) = logio(Mi) + /31ogio( A/star /A/o) 
1 (A4tar/Mo)* 1 
l + {A/Eitar/A/o)"^ 2 




Normalization 


logi(,(A/i) = 12.35 - 0.30(1 - a)* 


login (Afi) = 0.03502^ - 


0.1922 + 10.199 




Characteristic Mass 


logio(A/o) = 10.72 - 0.59(1 - a)* 


logio(A/2) = 0.005092^ -t 


- 0.002992 + 11.824 




Low Mass Slope 


P = 0.43 - 0.18(1 - a)* 


a = -0.20762^ +0.7522 


+ 2.423 




High Mass Slope 


5 = 0.56 - 0.18(1 - a)* 
7 = 1.54 - 2.52(1 - a)* 


/3 = 0.1202^ - 0.09942 + 


0.206 




Stewart et al. (2009b) 


Mgas(Mstar) ^ 


Mstar(MgE 






Fitting Function 


Af^as — 0.04Afstar ( ^ 


/ A/ 

A.r _ OK » f / ^'-'gas ' 




i«gas ^■"^"^^^^'^ \4,.5xlO'-'- Mq J 


x.-islar jvjgas 






Redshift dependence 


M = 0.59(1 + 2)0« 


M = 0.59(1 + 2)0-« 


Scatter 


(7(logio(Afgaa/Afstar)) = 0.34 - 0.191og(l + z) 


(T(logiQ(Afstar/Afgas)) = 


0.34 - 0.191og(l + z) 





t These fits are "original" fits given by the corresponding paper. 

•I- These fits are newly derived from the original data (or fitting functions). 

*a is the scale factor, with a = 1 being present day, and a = 1/(1 + z). 
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Abundance 
Matching and 
Gas Model 
Options 



Input 
Parameters 

(Empt/ field "hints" 
show value I'anges 
foi valid results! 



Switch Input 
Screens 



■Gdtjxy Mass Congr^r^lOTl 



Abundance Mdtchinc Model: 

i^^ Xe:!- -Id I 



lAllQW i(i 



^jMear value | Allow sc^niAr 



Lfig MvU 



H^dStilfl 



Log MWB = 0.0 
Log UgX, ^ 0.0 

Log MlMf = 0,0 
(j« Feat. = 0.0 



Output 
Field 



CilcuJait 


1 kMr tc... 


^ Mstar 10... 


Mgasto... t1c#r 



^^^^ 



site direction, priority should be given to Mstar(-^^vii ) for 
their model. 

Despite these small deviations, the errors introduced 
by fitting the various models to the functional forms pro- 
vided are no larger than the typical uncertainties inher- 
ent in the abunda.nce m atching technique 0.25 dex; 
see lBehroozi et al]|2010[) . While I refer the reader to BIO 
for a thorough exploration of the statistical uncertainties 
and errors associated with abundance matching models, 
I do caution that all of these models are inherently less 
certain for z > 1. due to uncertainties in the observed 
stellar mass function at higher redshifts. 

3. GALAXY GAS MODEL 

In order to estimate the gas content of galaxies as 
a function of stellar mass, I adopt the relation from 
IStewart etlD (j2009bf ). henceforth S09, which quanti- 
fies observationally-inferred relations between galaxy gas 
fraction and stellar ma ss. Specifically, S 09 characterizes 
gas fraction data from iMcGaughl 120051 (di sk-dominated 
galaxies at z = 0) and lErb et al. 020061 (UV-selcctcd 
galaxies at z 2) with a relatively simple function of 
stellar mass and redshift (their Equation 1). 

( 4.5 x ^y^Mj 

As with the abundance matching fits, this relation also 
depends on redshift. Though very high stellar mass 
galaxies tend to have similar gas content at z = — 2, 
the slope of this relation is much steeper at higher red- 
shift, a consequence of galaxies being more gas-rich at 
earlier times. The redshift evolution of /i(z) is given in 
the left column of Table |TJ Unlike abundance matching, 
the relation between stellar mass and gas mass is consid- 
erably more stochastic, with significant intrinsic scatter 
(especially at low redshift). While the above equation 
captures the mean relation, S09 also characterizes the 
typical scatter about log(Mstai7Afgas), characterized (in 
log space) by a gaussian with standard deviation a) . This 
log scatter is given in the left column of Table HJ 

Because the relation between Mgtai- and Afgas is signif- 
icantly less complex than that between A'/gtar and Mvir, 
it is a simple exercise to invert Equation [5] in order to ob- 
tain Mstar (Afgas)- Similarly, because the scatter is char- 
acterized by the log of the ratio Afstar/Afgas, the same 
equation for the magnitude of scatter is also applicable 
to the log of the inverse ratio, Afgas/Afstar (see the right 
column of Table |T]) . 

3.1. Gas Model Discussion 

S09 finds that the above characterization, which is de- 
fined by the gaseous propert ies of disk galax ies at z = 
and z = 2 (|McGaughl 120051 lErb et al.ll2006l ) is also re- 
mark ably consistent with other z = measurements 
(e.g.. lKannappanll2004l:lWei et al.ll2010D. as well as inter- 
mediate re dshift observations (e.g.. iWrig ht et al.l I200I 
lHamme~t al. 200^. GalMass provides two methods 
for utilizing this galaxy gas model. The first option is 
to only return the average gas mass for a disk galaxy 
of a given stellar mass. The second option, labeled "al- 
low scatter" uses a pseudo-random number generator to 
include scatter (up to ±ct) from the mean value. 



Figure 1. Screcnshot of the GalMass application. The upper 
(blue) region shows the abundance matching and galaxy gas model 
options (see ^and Q. The left (green) and right (purple) regions 
show the input and output fields, respectively. When the input 
boxes are empty, they display "hint" values that span typical mass 
and redshift ranges where the models are applicable. The "Calcu- 
late" button reads from the input boxes, applies the proper fitting 
functions, and updates the output field. The bottom (red) region 
allows the user to switch screens, changing the input mass parame- 
ter from between Afvir , A^star , and Mgas . The bottom right "Clear" 
button clears all screens of their input and output values. 

Note that this gas model is based on low redshift sam- 
ples of isolated disk galaxies. It does not claim to repro- 
duce realistic gas content for satellite galaxies that may 
have had gas stripped from their halos upon infall, nor 
is it designed to reproduce the gas content of elliptical 
galaxies. Specifically, the model is designed to fit ob- 
servations in the stellar mass range lO^'^Af© < Afgtar < 
10^^-^ Mq. While GalMass will allow for input parame- 
ters outside of this range (extrapolating the fitting func- 
tion to higher or lower masses), these resulting values 
will not be robust. 

As with the abundance matching data, this relation 
is also les s relia ble at higher redshift, as the data from 
lErb et all (|2006f ) is based on indirect estimates of galaxy 
gas fractions. The maximum redshift at which GalMass 
will produce outputs based on this gas model is z = 2. 

4. IMPLEMENTATION 

Throughout the application, masses are always given 
in units of solar mass, with log values displayed for com- 
pactness. Halo virial mass, galaxy stellar mass, galaxy 
gas mass, and galaxy baryonic mass (gas and stars) are 
denoted "Mvir," "Mstar," "Mgas," and "Mbar" respec- 
tively. The application also computes the galaxy gas frac- 
tion ("Gas Frac") defined as Afgas/(Afgas + A/star)- 

The basic layout of GalMass consists of three inter- 
changeable screens, depending on the input parameter 
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desired by the user. The buttons that switch between 
screens are found at the bottom of each screen. For ex- 
ample, the "Mvir to..." button corresponds to the screen 
where the halo virial mass is expected as an input, and 
the galaxy's stellar and gas masses arc given as outputs 
(see Figure [T]). 

For a more detailed discussion of what mass regimes 
and redshifts abundance matching results can be trusted, 
I refer the reader to the original papers from whence 
I have derived fitting functions: CW09, MIO and BIO. 
Still, to aid the user in knowing when the calculations 
are based on extrapolations beyond well-tested regimes, 
GalMass displays a warning message (but output fields 
are still populated) if any calculations result in one the 
following: Mvir > W^^Mq; Mvir < 10^°Mq; or M^tar < 
IO^Mq. a more severe warning message is displayed 
(and output fields are cleared to null values) if any input 
or output values contain the following: Mvir > lO^^M©; 
Mstar > 1O^^M0; z > 2; or any mass < IM0. 

Lastly, I re-emphasizc that both the abundance match- 
ing data and the galaxy gas models, as implemented here 
in GalMass, are designed for isolated galaxies. The gas 
model will almost certainly over-estimate the gas con- 
tent of satellite galaxies, and abundance matching mod- 
els typically match stellar masses to dark matter virial 
masses upon first infall (when the galaxy was last iso- 
lated). 
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